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Abstract 

The field of TeV gamma-ray astronomy has produced many exciting results over the last decade. Both the source 
catalogue, and the range of astrophysical questions which can be addressed, continue to expand. This article presents 
a topical review of the field, with a focus on the observational results of the imaging atmospheric Cherenkov telescope 
arrays. The results encompass pulsars and their nebulae, supernova remnants, gamma-ray binary systems, star forming 
regions and starburst and active galaxies. 
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1. Introduction 

Teraelectronvolt (TeV) astronomy concerns the study 
of astrophysical sources of gamma-ray photons, with 
energies in the range between ~ 30 GeV and ~ 30 TeV. 
The TeV range is one of the most recent windows of 
the electromagnetic spectrum to be opened for study, 
beginning with the identification of the first source, the 
Crab Nebula, in 1989 The results have been im- 
pressive, with significant advances in instrumentation 
leading to the detection of well over 100 sources over 
the last decade. The goals of TeV astronomy are wide- 
ranging, but can broadly be described as the study of 
sites of relativistic particle acceleration in the Universe, 
both hadronic and leptonic. This encompasses a huge 
range of size scales and energetics, from the interactions 
of galaxy clusters, to the magnetospheres of individual 
pulsars. 

Numerous recent reviews of TeV gamma-ray astron- 
omy have been written (e.g. 12, S 3]), but the field is 
rapidly evolving. For example, since the extensive re- 
view by Hinton & Hofmann in 2009 | 2], the number of 
sources has grown from around 80 to more than 120, 
and a number of new source classes have been identi- 
fied. This paper therefore aims to provide an update, 
and to supplement the existing reviews with a summary 
of the current observational status. The primary focus 
of the review is on the results from imaging atmospheric 
Cherenkov telescopes, with a relatively brief discussion 
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of the air shower particle detector experiments. We use 
the definitions outlined by Aharonian [5] as follows: 
"high energy", or "GeV", astronomy refers to the en- 
ergy range from 30 MeV to 30 GeV, while "very high 
energy", or "TeV", astronomy, refers to the range from 
30 GeV to 30 TeV. 



2. A Brief History 

The development of ground-based gamma-ray as- 
tronomy is closely linked to the study of cosmic rays 
and cosmic ray air showers. The idea of searching for 
astrophysical gamma-ray sources at ~ 100 MeV ener- 
gies was first proposed by Morrison in 1958 [6]. A pre- 
diction of a very high TeV gamma-ray flux from var- 
ious sources, including the Crab, which might be de- 
tectable with an air shower particle array at high alti- 
tude, was made by Cocconi in 1959 [7]. Cherenkov ra- 
diation associated with large cosmic ray air showers was 
first detected by Galbraith & Jelley in 1953 f?], and the 
possibility of using this phenomenon to study gamma- 
ray initiated showers led to the development of a num- 
ber of dedicated facilities in the 1960s. This eff'ort was 
boosted by the apparent detection of a gamma-ray signal 
from the black-hole binary Cygnus X-3 by both particle 
air shower arrays and atmospheric Cherenkov detectors. 
With hindsight, this detection was likely spurious, as 
were numerous unsubstantiated claims throughout the 
1970s, 1980s and into the 1990s. 

The field reached a firm experimental footing with the 
development of the imaging technique, which provides 
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a method of effectively discriminating between gamma- 
ray initiated showers and the background of cosmic ray 
showers, based on the morphology of their Cherenkov 
images, and guided by the results of Monte Carlo simu- 
lations f5,To]. This technique was applied by the Whip- 
ple collaboration to detect steady gamma-ray emission 
from the Crab Nebula using a 10m reflector with a 37- 
element photomultiplier tube camera in 1989 idl]. A 
number of imaging atmospheric Cherenkov telescopes 
(lACTs) were subsequently developed around the world 
(including Durham, CANGAROO, Telescope Array, 
Crimean Astrophysical Observatory, SHALON, TAC- 
TIC), with the northern hemisphere instruments (Whip- 
ple, HEGRA and CAT) leading the field. The 1990s saw 
two particularly important developments: the detection 
of the first extragalactic sources by the Whipple Collab- 
oration, starting with the nearby blazars Markarian 421 
fill] and Markarian 501 1.121 , and the application of the 
stereo imaging technique by the HEGRA array lITsll . 
HEGRA consisted of 5 telescopes of modest aperture 
(< 10 m^), and demonstrated that the combination of 
Cherenkov image information from multiple telescopes 
located within the same Cherenkov light pool could dra- 
matically improve the sensitivity of the technique. 

Despite this progress, the relative scarcity of bright 
TeV gamma-ray sources (< 10 were identified by 2000) 
highlighted the necessity for improved instrumentation. 
Cherenkov wavefront samplers such as CELESTE and 
STACEE attempted to probe to lower energies, and 
hence higher gamma-ray fluxes and larger distances, 
using converted solar farms; however, the difficulty of 
discriminating gamma-rays from the cosmic ray back- 
ground using this technique limited its effectiveness. 
Successful gamma-ray observations using a particle de- 
tector were made by the Milagro experiment, which ran 
from 2000 to 2008, providing a survey of the north- 
ern sky with modest sensitivity. Starting with the com- 
missioning of H.E.S.S. in 2003, the new generation of 
lACTs - H.E.S.S., MAGIC and VERITAS - have pro- 
vided the required order of magnitude improvement in 
sensitivity, and firmly established gamma-ray studies as 
an important astronomical discipline. 

3. Current Status and Instrumentation 

3.1. Imaging Atmospheric Cherenkov Telescopes 

The Cherenkov emission from an air shower forms a 
column of blue Ught in the sky, with the maximum emis- 
sion occuring around 10 km above sea level at TeV ener- 
gies. Cherenkov telescopes used to record these images 
are essentially rather simple devices, consisting of a 



large, segmented optical flux collector used to focus the 
Cherenkov light onto an array of fast photo-detectors. 
The optical specifications are not terribly strict; an op- 
tical point spread function width of typically ~ 0.05° 
is adequate. The photo-detector array (usually < 1000 
photomultiplier tubes) comprises a crude camera, cover- 
ing a few degrees on the sky. A Cherenkov flash triggers 
read-out of the photo-detectors, with a read-out window 
defined by the timescale of the arrival of the Cherenkov 
photons (~ 10 ns). 

There are currently three major imaging atmospheric 
Cherenkov telescope systems in operation. H.E.S.S., 
located in the Khomas Highland of Namibia (-23°N, 
-16°W, altitude 1800 m), consists of four telescopes ar- 
ranged on a square with 120 m side length. Each tele- 
scope has a mirror area of 107 m^ and is equipped with 
a 960 pixel camera covering a 5° field of view. VER- 
ITAS, at the Fred Lawrence Whipple Observatory in 
southern Arizona (32°N, 111°W, altitude 1275 m) has 
similar characteristics, with 4 telescopes of 107 m^ area 
and 499-pixel cameras, covering 3.5°. MAGIC (28°N, 
17° W, altitude 2225 m) originally consisted of a single, 
very large reflector (236 m^) on the Canary island of La 
Palma, with a 3.5° camera. In 2009, a second telescope 
with the same mirror area was installed at a distance of 
85 m from the first. 

Each of these facilities work in a similar fashion. 
Cherenkov images of air showers are recorded at a rate 
of a few hundred Hz, and analyzed offline. The over- 
whelming majority of these images are due to cosmic 
ray initiated air showers. Gamma-ray showers can be 
discriminated from this background based on the image 
shape and orientation. Gamma-ray images result from 
purely electromagnetic cascades and appear as narrow, 
elongated ellipses in the camera plane. The long axis of 
the ellipse corresponds to the vertical extension of the 
air shower, and points back towards the source position 
in the field of view. If multiple telescopes are used to 
view the same shower, the source position is simply the 
intersection point of the various image axes (illustrated 
schematically in Figure [T]l. Cosmic -ray primaries pro- 
duce secondaries with large transverse momenta, which 
initiate sub-showers. The Cherenkov images of cosmic- 
ray initiated air showers are consequently wider then 
those with y-ray primaries, and form an irregular shape, 
as opposed to a smooth ellipse. In addition, since the 
original charged particle has been deflected by galactic 
magnetic fields before reaching the Earth, cosmic-ray 
images have no preferred orientation in the camera. 

Cherenkov light reaching the ground from air show- 
ers peaks at optical/UV wavelengths, and so lACTs op- 
erate only under clear, dark skies. Both MAGIC and 



2 



3.2. Particle Detectors 
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Figure 1 : A schematic illustration of an atmospheric Cherenkov tele- 
scope array. The primary paiticle initiates an air shower, resulting in a 
cone of Cherenkov radiation. Telescopes within the Cherenkov light 
pool record elliptical images; the intersection of the long axes of these 
images indicates the arrival direction of the primary, and hence the 
location of a y-ray source in the sky. 



VERITAS have demonstrated that useful observations 
can be made when the moon is visible above the hori- 
zon, but the typical duty cycle of these instruments is 
still limited to ~ 1200 hours per year (< 15%). Given 
the small field of view of lACTs, regions of the sky 
containing one or more source candidates are usually 
targeted for observations. Surveys can only be accom- 
plished slowly, by tiling regions of the sky with overlap- 
ping fields-of-view. The sensitivity of the current gen- 
eration of lACTs is sufficient to detect the Crab Nebula 
in under a minute, and a source with 1% of the Crab 
flux (~ 2 X 10"'^ m-^ s^' above 1 TeV) in ~ 25 hours. 
The angular and energy resolution of the technique are 
energy-dependent, with typical values of < 0.1° and 
< 15% per photon, respectively, at 1 TeV. 

The catalog of TeV sources grew rapidly with the 
commissioning of H.E.S.S. in the southern hemisphere, 
which provided the first high sensitivity observations of 
the densely populated inner Galaxy. It has continued to 
expand in recent years as MAGIC and VERITAS have 
come online, and now numbers 130 sources, as listed in 
the online catalog TeVCat 1 14] 0. Figure |2] shows the 
locations of these sources in Galactic coordinates. 



'We note that a small fraction of these 130 sources may be dupli- 
cate detections of the same object, while other detections remain con- 
tentious or unconfirmed. Full details are available in the individual 
source annotations in TeVCat, available at http://tevcat.uchicago.edu 



The direct detection of air shower particles using ar- 
rays of particle detectors at ground level offers some 
important advantages over the atmospheric Cherenkov 
technique. In particular, observations can be made con- 
tinuously, both day and night, and over the entire view- 
able sky (a field-of-view of > 1 steradian). These ad- 
vantages are off'set, however, by the rather low sensitiv- 
ity to point sources, which is primarily due to the dif- 
ficulty of rejecting the substantial background of cos- 
mic ray initiated air showers. The angular and en- 
ergy resolution of these detectors are also significantly 
worse than lACTs. Early claims for emission from bi- 
nary systems using sparse particle detector arrays were 
not confirmed by later, more sensitive instruments (e.g. 
lll5 . 16]), indicating the need for a new approach to this 
problem. Milagro, which operated between 2000 and 
2008 in northern New Mexico (36°N, 107°W), was the 
first successful attempt at this. 

The Milagro detector consisted of a large water reser- 
voir (60 X 80 m) at an altitude of 2630 m, covered with 
a light-tight barrier, and instrumented with PMTs. The 
central reservoir provided high-resolution sampling of 
air shower particles over a relatively small area (com- 
pared to the air shower footprint). In 2004 an array 
of 175 small outrigger tanks were added, imgularly 
spread over an area of 200 x 200 m around the cen- 
tral reservoir This configuration, coupled with the de- 
velopment of analysis techniques for cosmic ray back- 
ground discrimination, provided sufficient sensitivity 
for the first comprehensive survey of the northern TeV 
sky. The results showed strong detections of the bright, 
known TeV sources Markarian 421 and the Crab Neb- 
ula, along with the detection of three extended sources 
in the Galactic plane, each with integrated fluxes com- 
parable to the Crab Nebula at 20 TeV. A few less sig- 
nificant source candidates were also identified in the 
plane, and a reanalysis following the launch of Fermi- 
LAT also showed fourteen 3cr excesses co-located with 
bright Galactic LAT sources [17]. The Milagro results 
for the region around the Galactic plane are shown in 
Figure [3] 

Milagro ceased operation in 2008; however, two large 
particle detector arrays remain in operation at very high 
altitude in Tibet. ARGO-YBJ is located at 4300 m 
in Yangbajing, and consists of a single layer of re- 
sistive plate chambers completely covering an area of 
110 X 100 m. The results of 1265 days of observa- 
tions were recently presented fl^, showing > 5cr detec- 
tions of the Crab Nebula, Markarian 421 and two Mila- 
gro sources. The fact that one of the brightest Milagro 
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Figure 3: The Milagro survey of the Galactic plane. The z-axis is the 
pre-trials statistical significance, with a fixed maximum of Icr. Figure 
from Abdo et al. \\^. 



sources, MGRO J2019+37, is not detected in these ob- 
servations presents something of a mystery. The Tibet 
ASy air shower array, also at Yangbajing, consists of 
~ 750 closely-spaced scintillation detectors covering an 
area of 36900 m^, and has demonstrated that this tech- 
nique is also practical for the detection of bright TeV 
sources |20]. 



4. Extragalactic TeV Sources 

4.1. Blazars 

Approximately 1 % of all galaxies host an active nu- 
cleus; a central compact region with much higher than 
normal luminosity. Around 10% of these Active Galac- 
tic Nuclei (AGN) exhibit relativistic jets, powered by 
accretion onto a supermassive black hole. Many of the 
observational characteristics of AGN can be attributed 
to the geometry of the system; in particular, the orien- 
tation of the jets with respect to the observer Blazars, 
which host a jet oriented at an acute angle to the line 
of sight, are of particular interest for gamma-ray astron- 
omy, as the emission from these objects is dominated by 
relativistic beaming effects, which dramatically boost 
the observed photon energies and luminosity. 

The first extragalactic source discovered at TeV ener- 
gies was Markarian 421 1 1 1], a blazar of the BL Lacer- 
tae sub-class. The extragalactic TeV catalog now com- 
prises ~ 50 objects, and continues to increase steadily 
(Figure nil. Blazar SEDs show a double-peaked struc- 
ture in a vFy representation of their spectral energy dis- 
tribution (SED), with the lower frequency peak usu- 
ally attributed to synchrotron emission of energetic elec- 
trons, and the higher frequency peak to inverse Comp- 
ton. BL Lac objects are further classified as low-, 
intermediate- or high-frequency peaked, according to 
the location of the peak of their synchrotron emission. 
The majority (~ 80%) of the known TeV blazars are 
high-frequency peaked objects, in part because of inher- 
ent biases in the target selection: initially, objects were 
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Figure 4: The location of the 49 AGN (BL Lac, FSRQs and radio galaxies) with detected TeV emission (from TeVCat [IJ]). On the left is the 
spatial distribution, in Galactic coordinates; The right plot shows the redshift distribution, for 41 of the objects with redshifts listed in the literature. 



chosen based primarily upon their radio and X-ray spec- 
tral properties (e.g. nllv ). More recently, the TeV obser- 
vatories have expanded their selection criteria, using ad- 
ditional guidance in the form of FermZ-LAT results, and 
multi-wavelength observation triggers. This has broad- 
ened the catalog to include examples of intermediate- 
and low-frequency peaked objects. The overall data 
quality has also improved markedly since the launch of 
Fermi: Figure |5] shows a recent compilation of spectral 
measurements for the bright TeV blazar, Markarian 501 
1 2211 . taken during an extensive multi wave length cam- 
paign in 2009. 

The mechanisms which drive the high energy emis- 
sion from blazars remain poorly understood, and a full 
discussion is beyond the scope of this review. Briefly; in 
leptonic scenarios, a population of electrons is acceler- 
ated to TeV energies, typically through Fermi accelera- 
tion by shocks in the AGN jet. These electrons then cool 
by radiating X-ray synchrotron photons. TeV emission 
results from inverse Compton interactions of the elec- 
trons with either their self-generated synchrotron pho- 
tons, or an external photon field. The strong correla- 
tion between X-ray and TeV emission which is often 
observed provides evidence for a common origin such 
as this, although counter-examples do exist ll23ll . An- 
other class of models has hadrons as the primary particle 
population, which can then produce TeV gamma-rays 
through subsequent interactions with target matter or 
photon fields. Hadronic models are less favoured, typi- 
cally, in part because the cooling times for the relevant 
processes are long, making rapid variability difficult to 
explain. One exception to this is the case of proton syn- 
chrotron emission, which may provide a plausible al- 
ternative, in which the emission results from extremely 
high energy protons in highly^magnetized (B ~ 100 G), 
compact regions of the jet |24]. 

Many of the AGN detected at TeV energies exhibit 
extreme variability. The timescales can range from 
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Figure 6: Integrated flux (> 200 GeV) versus time for H.E.S.S. ob- 
servations of PKS 2155-304 on MJD 53944 (28 July, 2006). The data 
are binned in 1 -minute intervals, and the horizontal dashed line shows 
the steady flux from the Crab Nebula for comparison. Figure from 



years to minutes, and the observed flux can change by 
more than an order of magnitude. Figure |6] shows the 
H.E.S.S. lightcurve from July 2006 for one of the most 
extreme examples, the BL Lac object PKS 2155-304 
nisi 12511 . Such rapid variability can be used to place 
constraints on the size of the emission region, which de- 
pend upon the Doppler factor, S. 6 itself is constrained 
by the requirement that the emission region should be 
transparent to gamma-rays (e.g. [26]). Extremely rapid 
TeV gamma-ray variability of distant blazars can also 
be used to place limits on the energy-dependent viola- 
tion of Lorentz invariance [27, 28], as predicted in some 
models of quantum gravity. 

Related to the BL Lacertae objects are Flat Spectrum 
Radio Quasars (FSRQs). These are characterized pri- 
marily by their intense UV emission, associated with an 
accretion disk, strong broad emission lines in the opti- 
cal spectrum, and infra-red emission associated with a 
dusty torus. FSRQs are similar to low-frequency peaked 
BL Lacs, in that the X-ray emission is dominated by 
the inverse Compton peak of the SED. Despite this, in- 
verse Compton emission can extend up to TeV ener- 
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Figure 5: Extensive multifrequency measurements showing the spectral energy distribution of Markarian 501 for observations in 2009. Emission 
from the host galaxy is clearly visible at infrared/optical frequencies. The VERITAS data are divided to show both the average spectrum (red 
circles), and the spectrum during a 3-day flare (green tiiangles). Figure from Abdo et al. t22ll 



gies, particularly during intense flaring episodes, and 
three FSRQs have recently been detected by lACTs 
fPKS 1222+21 ii , PKS 1510-089 ^ and 3C279 

iH). 



A unique additional case is the TeV detection of 
IC 310, a "head-tail" radio galaxy in the Perseus clus- 
ter, possibly hosting a low-luminosity BL Lac nucleus 
1 3311 . Head-tail radio galaxies display a distinctive radio 
morphology, consisting of a bright "head" and a fainter 
"tail", which is believed to be the result of their rapid 
motion with respect to the intracluster medium. The 
source was originally identified as a VHE emitter in an 
analysis of the highest energy (> 30 GeV) Fermi pho- 
tons by Neronov et al. fsT], and then subsequently de- 
tected from the ground by MAGIC [,35]. Neronov et al. 
considered the intriguing possibility that the emission 
might originate with particles accelerated at the bow 
shock formed by the interaction between the relativistic 
outflow from the galaxy and the intracluster medium. 
This scenario is ruled out by the detection of variability 
in the TeV flux by MAGIC, and the more familiar BL 
Lac mechanisms are now favoured. 

Figure |4] also shows the distribution of measured red- 
shifts for 41 TeV AGN. Some of these measurements 
are rather uncertain, since BL Lac optical spectra, by 
definition, do not contain strong emission lines. The 



most distant object detected is 3C279 fy^, with a rel- 
atively modest redshift of z = 0.5362. The population 
is truncated at large distances due to the absorption of 
TeV gamma-rays by electron-positron pair production 
with the low energy photons of the extragalactic back- 
ground light (EBL). This effect is energy dependent, and 
can strongly modify the observed VHE spectra of extra- 
galactic sources. While this limits the observation of 
more distant TeV sources, it also provides a mechanism 
by which to infer the intensity of the EBL, using rea- 
sonable assumptions about the intrinsic TeV spectra at 
the source ll36ll . The EBL provides a calorimetric mea- 
sure of the complete history of star and galaxy forma- 
tion in the Universe, but is difficult to measure directly, 
due to the presence of bright local foreground sources 
of emission, such as zodiacal light. Presently, all of 
the TeV blazar measurements are consistent with a rel- 
atively low level of EBL, with the constraints derived 
from VHE measurements now approaching the lower 
limits derived from galaxy counts ll37ll38ll39ll . 



TeV blazar observations have also been suggested as 
probes of other physical phenomena, such as the ac- 
celeration and propagation of ultra-high energy cosmic 
rays 14^ 41], or, more speculatively, the production of 
axion-like particles |42, 43]. Various authors have also 
discussed the possibility that TeV observations may be 
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used to measure or constrain the strengfli of the inter- 
galactic magnetic field (IGMF) (e.g. H Si)- Tem- 
poral, spectral and spatial signatures of the IGMF are 
all possible; however, the fact that blazars are intrin- 
sically variable gamma-ray sources limits the power of 
this technique. Accounting for this, Dermer et al. derive 
a lower limit of B,cmf > 10"* G H. 



4.2. Radio Galaxies 

As described above, the TeV fluxes from blazars are 
dramatically enhanced by the effects of Doppler boost- 
ing. Nearby radio galaxies, in which the jet is not di- 
rectly oriented towards the line-of-sight, provide an al- 
ternative method by which to investigate the particle 
acceleration and gamma-ray emission from relativistic 
outflows in AGN. The advantage of studying such ob- 
jects lies in the fact that the jets can be resolved from 
radio to X-ray wavelengths, allowing the possibility of 
correlating the gamma-ray emission state with observed 
changes in the jet structure. Three radio galaxies have 
been identified as TeV emitters: M 87, Centaurus A and 
NGC 1275. 

M 87 is the most well studied of these, and was first 
reported as a gamma-ray source by the HEGRA collab- 
oration ('46*1, with subsequent confirmation by H.E.S.S. 

VERITAS 1 48] and MAGIC L49J. M 87 is a giant 
radio galaxy at a distance of 16.7 + 0.2 Mpc, displaying 
a prominent misaligned jet, with an orientation angle of 
< 20° to the line-of-sight. The mass of the central black 
hole is estimated to be ~ 3 x 10^ Mq. The TeV source 
is strongly variable, and has undergone three episodes 
of enhanced emission in 2005, 2008 and 2010 (Fig- 
ureEJ. The results are summarized by Abramoswski et 
al ll50l] . Causality arguments use the shortest variability 
timescale of around 1 day to place strong constraints on 
the size of the TeV emission region, corresponding to 
only a few Schwarzschild radii. The TeV emission re- 
gion cannot be dkectly resolved with lACTs, but corre- 
lations with contemporaneous X-ray and radio observa- 
tions provide some clues to its location. Two structures 
are of particular interest: the inner region close to the 
central black hole (the "core"), and HST-1, a bright jet 
feature first resolved in the optical band by the Hubble 
Space Telescope. HST-1 underwent a multi-year flare in 
radio, optical and X-rays, peaking around the time when 
the first short-term variability was detected at VHE en- 
ergies Ii47.1 . In contrast to this, the 2008 VHE flare was 
accompanied by enhanced radio and X-ray fluxes from 
the core region fs^. The 2010 VHE flare showed no 
enhanced radio emission from the core, although an en- 
hanced X-ray flux was observed 3 days after the VHE 
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Figure 7: VHE light curves of the three flares from M87 observed 
in 2005, 2008 and 2010, showing integral fluxes above an energy of 
350 GeV. Figure from |50]. 



peak. The results, therefore, remain somewhat am- 
biguous, and the possibility remains that the observed 
VHE flares may have different origins, or that the to- 
tal VHE emission may be the sum of multiple compo- 
nents. Given the burdensome observing requirements 
for instruments with a limited duty cycle, M 87 repre- 
sents the best example of the importance of data-sharing 
and coordinated observing planning between the vari- 
ous lACTs. 

The closest active galaxy, Centaurus A, at a distance 
of 3.8 Mpc, was identified as a VHE source in a deep, 
120 hour exposure by H.E.S.S. |52]. Cen A is among 
the faintest VHE sources detected, with a flux of 0.8% 
of the Crab Nebula above 250 GeV. The emission 
is steady, although the variability is not strongly con- 
strained as a result of the low flux level. As with M 87, 
numerous sites for the production of the TeV emission 
have been suggested, from the immediate vicinity of 
the central black hole (with a mass of ~ 5 x 10^ M©), 
to the AGN jet, or even beyond ['53]. The final, and 
most recent, addition to the known TeV radio galaxies 
is NGC 1275, identified by MAGIC Q. NGC 1275 
is the central galaxy of the Perseus cluster, at a distance 
of 72.2 Mpc. The VHE emission has a flux of ~ 3% 
of the Crab Nebula above 100 GeV and exhibits a very 
soft spectrum, with a power-law index of 3.96 + 0.37 
1.55,1 . indicating a sharp turnover from the measured 
Fermi spectrum at lower energies, where the index is 
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2.09 + 0.02 115611 ■ 



4.3. Starburst Galaxies 



Starburst galaxies are those which exhibit an ex- 
tremely high rate of star formation, sometimes triggered 
by interaction with another galaxy. High cosmic -ray 
and gas densities in the starburst region make these 
objects promising targets for gamma-ray observations, 
with emission predicted to result from the interactions 
of hadronic cosmic rays in the dense gas. TeV emis- 
sion has now been identified from two starburst galax- 
ies: M 82 [57] and NGC 253 [58]. 

M 82 is a bright galaxy located at a distance of ap- 
proximately 3.9 Mpc, with an active starburst region at 
its centre. The star formation rate in this region is ap- 
proximately 10 times that of the Milky Way, with an 
estimated supernova rate of 0. 1 to 0.3 per year A deep 
VERITAS exposure (137 hours) in 2008-2009 resulted 
in a detection of gamma-ray emission from M 82 with 
a flux of (3.7 + 0.8,,fl, + 0.7„,„) x 10"'^ cm^^ s"' above 
the 700 GeV energy threshold of the analysis. NGC 253 
lies at a distance of 2.9 - 3.6 Mpc, and also has a cen- 
tral, compact (~ 100 pc) starburst region. The super- 
nova rate in this region is estimated at ~ 0.03 per year. 
TeV emission was detected by H.E.S.S. with an inte- 
grated flux above 220 GeV of (5.5 + l.Qstat ± ^-^sxst) x 
IQ-^^ cm-2 m. 

The emission from both M 82 and NGC 253 is con- 
sistent with the predictions of models based on the ac- 
celeration and propagation of cosmic rays in the star- 
burst core (e.g. ll59ll ). assuming that they act as efficient 
"proton calorimeters" (i.e. cosmic rays lose the major- 
ity of their energy to collisions). In this case, the esti- 
mated cosmic ray density in the starburst region is 2-3 
orders of magnitude larger than that of the Milky Way. 
An alternative explanation is proposed by Mannheim et 
al. ll60ll . who suggest that the TeV luminosity is consis- 
tent with the combined emission from a large population 
of pulsar wind nebulae, which result from the elevated 
supernova rate. More accurate TeV spectra, and obser- 
vation of other starburst classes, such as ultra-luminous 
infrared galaxies, should provide more insight in the fu- 
ture. 



4.4. The Large Magellenic Cloud 

Galaxies of the Local Group are also of interest to 
TeV observatories, although the predicted fluxes due to 
cosmic ray acceleration and propagation lie below the 
current instrumental sensitivity in the TeV range. Given 



their proximity, the possibility arises of detectable emis- 
sion from individual objects, or from localized re- 
gions, particularly in the Milky Way's satellite galax- 
ies. H.E.S.S. has recently identified the first such ob- 
ject in the Large Magellenic Cloud (LMC), at a dis- 
tance of 48 kpc [61] . An unresolved source with a flux 
of ~ 2% of flie Crab Nebula (1.5 x lO^'^ erg cm^^ s"' 
between 1 and 10 TeV) was detected, consistent with 
the location of PSR J0357-6910. This object is the 
most powerful pulsar known, with a spin-down energy 
of E= 4.8 X lO"*^ erg s"'. Given the positional coinci- 
dence, and based on comparisons with similar objects 
within our Galaxy, it seems Ukely that the TeV emis- 
sion is due to inverse Compton emission from electrons 
in the pulsar's wind nebula interacting with a strong in- 
frared target photon field. This is in contrast to the ex- 
tended GeV emission which has been observed from the 
LMC by Fermi-LAI, which is attributed to cosmic ray 
acceleration and interactions in the massive star forming 
region of 30 Doradus 1162 1. 



4.5. Other Extragalactic TeV targets 

TeV observations have also been used to place im- 
portant constraints on the gamma-ray emission from nu- 
merous undetected extragalactic source classes, includ- 
ing galaxy clusters and potential sources of ultra-high 
energy cosmic rays. Here we summarize two of the 
most important: gamma-ray bursts, and the predicted 
sites of dense regions of dark matter particles. 
Gamma-Ray Bursts: GRBs are the signatures of brief, 
extremely energetic explosions which occur at cosmo- 
logical distances. They are observationally divided into 
short and long classes, which are presumably the result 
of diff'erent progenitor systems. Long duration GRBs 
are generally ascribed to the collapse of massive, rapidly 
rotating stars into black holes. The origin of the short 
bursts is less certain, although neutron star - neutron 
star merger events are among the favoured candidates. 
The search for > 100 GeV emission has been a long- 
running goal of both the lACTs and particle detectors 
(see [63] for a review). Given the brief duration of GRB 
emission, the wide field-of-view of particle detectors is 
particularly important in this regard, although observa- 
tions so far have been hampered by limited sensitiv- 
ity. lACTs, conversely, must be re-pointed rapidly on 
receipt of an alert. The detection of delayed high en- 
ergy emission by Ferm/-LAT , lasting hundreds to thou- 
sands of seconds longer than the sub-MeV emission, 
has provided additional impetus to the search 11641 16511 . 
and lACTs now regularly target burst locations within 
< 100 s of the burst alert. The task is difficult, since 
the burst must also be at a small enough redshift such 
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Figure 8: Upper limits at the 95% confidence level on the velocity- 
weighted annihilation neutralino cross-section for different annihila- 
tion channels, based on VERITAS observations of the Segue I dwaif 
spheroidal galaxy. The dark band represents the typical range of pre- 
dictions. Figure from I73ll 



that the high energy emission is not completely sup- 
pressed by photon-photon pair production (z < 0.5). 
No convincing signals have been detected as yet (e.g. 
1 661, 67 , 68 , i69[l ). although predictions based on the 
brightest bursts observed by the LAT indicate that the 
potential for detecting TeV emission associated with a 
GRB is promising, assuming no intrinsic spectral cut- 
off of the high energy emission ItoIi . 



Extragalactic Dark Matter: The search for the self- 
annihilation signature of dark matter particles in astro- 
physical objects is wide-ranging, and complementary 
to direct detection techniques on the Earth (see 117 ill 
for an excellent review). As discussed later in this re- 
view, the centre of our own Galaxy is a natural tar- 
get for dark matter searches, and provides the most 
stringent limits to date ll72ll . Objects outside of our 
own Galaxy are also worthy of investigation, however, 
and are potentially much less affected by contamination 
from unknown astrophysical background sources (su- 
pernova remnants, pulsar wind nebulae, etc.). Dwarf 
spheroidal galaxies of the local group are among the 
most promising of these, due to their proximity and their 
presumed large dark matter content. The Sloane Digital 
Sky Survey has more than doubled the known popula- 
tion of dwarf spheroidals, providing additional targets 
for the TeV searches. No signals have been detected 
as vet, despite deep exposures on a number of objects 
III, 74, TSl 3 231 • Figure [8] shows upper limits on 
the annihilation cross-section derived for various anni- 
hilation channels using VERITAS observations of the 
Segue I dwarf spheroidal. 



5. Galactic TeV Sources 

There are presently ~ 80 known TeV sources within 
our Galaxy, as indicated either by their association with 
known Galactic sources at other wavelengths, or by 
their location in the Galactic plane - a particularly com- 
pelling argument when coupled with a resolvable angu- 
lar extent. The majority of these sources were identified 
as TeV emitters during the H.E.S.S. survey of the inner 
Galaxy 111- H.E.S.S. is the only lACT currently op- 
erating in the southern hemisphere, which allows it to 
view the inner Galaxy at high elevation, and hence with 
a low energy threshold and good sensitivity. H.E.S.S. 
was the first instrument with sufficient sensitivity to ob- 
serve sources with < 10% of the Crab Nebula flux in 
this region, and the results have been revelatory - the 
Galactic plane is densely populated with TeV sources, 
primarily clustered within the inner ±60° in Galactic 
longitude (Figure |9j. The most recent survey results 
consist of 2300 hours of observations, allowing the de- 
tection of over 50 sources within the range I - 280° to 
60° and h = -3.5° to -h3.5° \B\- Observations of the 
outer Galactic regions by VERITAS, MAGIC, Milagro 
and ARGO-YBJ have revealed a less densely populated 
sky, but containing some unique objects of particular in- 
terest for TeV studies. Many Galactic TeV sources are 
extended, allowing detailed studies of source morphol- 
ogy and spatially resolved spectra, while others are time 
variable and/or periodic. The various source classes are 
discussed in some detail below. 

5.7. The Galactic Centre and Rid^e 

A TeV source at the location of the Galactic Centre 



has been reported by various lACTs jsol Isil |82l [sl. 
Determining the nature of this source is a difficult task, 
due to the complexity of the region, which includes 
multiple different potential counterparts. The most de- 
tailed studies have been performed by H.E.S.S., which 
reveal that the emission is dominated by a bright cen- 
tral source, HESS J1745-290, lying close to the cen- 
tral supermassive black hole, Sgr A*. An additional, 
fainter, component is also seen, which extends in both 
directions along the Galactic plane llsill . The extended 
component is spatially correlated with a complex of gi- 
ant molecular clouds in the central 200 pc of the Milky 
Way, and the TeV emission can be attributed to the 
decay of neutral pions produced in the interactions of 
hadronic cosmic rays with material in the clouds. The 
central source is point-like, steady and exhibits a curved 
power-law spectrum fss'l. Its location with respect to 
three of the most likely counterparts is shown in Fig- 
ure do] This study reveals that the source centroid is 
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Figure 9: Significance map of the Galactic Plane from the original 
H.E.S.S. survey in 2004, based on 230 hours of observations. The 
complete survey now extends over the range from / = 280° to 60° and 
b = —3.5° to +3.5°, and comprises 2300 hours of observations f79ll . 
Figure from i78l1 



displaced from the radio centroid of the supernova rem- 
nant Sgr A East, excluding this object with high prob- 
ability as the dominant source of the VHE gamma-ray 
emission, and leaving Sgr A*and the pulsar wind nebula 
G359.95-0.04 as the most likely counterparts fse"]. 

The Galactic Centre is also a prime candidate region 
in which to search for gamma-ray emission due to dark 
matter particle self-annihilation. The analysis is com- 
plicated, however, because of the high background due 
to astrophysical sources. An analysis by H.E.S.S. using 
an optimized background subtraction technique shows 
no hint of a residual dark matter gamma-ray flux at a 
projected distance of r ~ 45 - 150 pc from the Galactic 
Centre ||72| . 



5.2. The Crab Nebula and Pulsar 

The Crab is the nearby (2.0 + 0.2 kpc) remnant of a 
historical supernova explosion, observed in 1054 A.D. 
There is no detected shell, and the broadband emission 
below ~ 100 MeV is dominated by a bright synchrotron 
nebula, powered by a central pulsar (PSR B0531-I-21). 
PSR B0531+21 is the most energetic pulsar in our 
Galaxy, with a pulse period of 33 ms, and a spin-down 
power of 4.6 x 10^*^ erg s The Crab Nebula and Pul- 
sar hold a unique place in the development of TeV as- 
tronomy: the birth of the field as an astronomical disci- 
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Figure 10: 90 cm VLA radio flux density map of the innermost 20 pc 
of the Galactic Centre, showing emission from the SNR Sgr A East. 
Black contours denote radio flux levels; the centre of the SNR is 
marked by the white square, and the positions of Sgr A* and G359.95- 
0.04 are given by the cross hairs and the black triangle, respectively. 
The 68% CL total error contour of the best-fit centroid position of 
HESS J1745-290 is given by the white circle. Figure from 186]: see 
that paper for full details. 



pline can be traced to the detection of the Crab Neb- 
ula TeV source by Weekes et al. using the Whipple 
10 m telescope, in the first application of the imaging at- 
mospheric Cherenkov technique ijjj]. Subsequently, the 
Crab has acted as a bright, standard candle for TeV ob- 
servatories. 

The SED of the non-thermal nebula emission dis- 
plays two components (Figure [TTJ. The dominant, low 
frequency component is explained by synchrotron ra- 
diation of high energy electrons spiraling in the mag- 
netic field of the nebula l87l 18811 . The higher frequency 
component is attributed to inverse Compton scattering 
of lower energy photons by these electrons, including 
microwave background photons, far infrared and the 
electron-synchrotron photons themselves. The electron 
population reaches energies of at least 10'^ eV, through 
acceleration occuring both in a relativistic particle out- 
flow driven by the spin-down energy of the pulsar, and 
in shocks where this outflow encounters the surround- 
ing nebula. The highest energy particles likely require 
an alternative explanation for their origin, such as direct 
acceleration in intense electric fields associated with the 
pulsar itself fso']. Observations of the synchrotron neb- 
ula from radio to X-ray wavelengths provide high reso- 
lution imaging of the emission region; however, the syn- 
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chrotron data alone only contain information concern- 
ing the product of the magnetic field strength and the 
relativistic electron density. Since the inverse Comp- 
ton component is independent of the magnetic field 
strength, the combined SED allows an estimate of the 
Nebula magnetic field, which is now constrained to be 
between 100 and 200 fiG H- 

The search for a VHE component to the pulsed emis- 
sion from the Crab has been long and, until recently, 
fruitless. Despite discouraging model predictions, and 
the detection of spectral cut-offs below 10 GeV in other 
pulsars, the fact that no super-exponential cut-off was 
observed by EGRET in the Crab Pulsar spectrum ini- 
tially provided some encouragement for a continued 
search by lACTs [91J. Fen«/-LAT subsequently ex- 
tended the GeV s pec trum and measured a sharp spectral 
cut-off at 6 GeV [90]. A campaign by MAGIC, using a 
specially designed "analog sum" hardware trigger, pro- 
vided the first ground-based measurement of gamma- 
ray emission from the Crab pulsar [j^. The initial 
MAGIC flux measurement above 25 GeV was, like the 
Fermi-LAI result, consistent with an exponential cut- 
off. The existence of an exponential cut-off is a natu- 
ral consequence of emission due to curvature radiation, 
as favored by various models (e.g. [|93|). Both VERI- 
TAS and MAGIC recently presented new results, which 
challenge this paradigm |94, 95, 96]. Pulsed emission 
is observed to extend up to well beyond 100 GeV, and 
the combined LAT-IACT spectrum can best be fit with a 
broken power law. The explanation for this high energy 
component is an open question, at present. Gamma-ray 
opacity arguments require that the emission zone of the 
highest energy photons must be at least 10 stellar radii 
from the surface of the neutron star - much further than 
previously assumed. The results require either a sub- 
stantial revision of existing models of high energy pul- 
sar emission, or the addition of a new component, not 
directly related to the MeV-GeV emission. 

As mentioned above, the Crab has been used as a 
standard candle in TeV astronomy, on the assumption 
that its emission was steady. This is now demonstra- 
bly false, at least at energies below ~ 1 GeV, with the 
detection of multiple day-scale flaring events |98l 89ll . 
and long term variation in the hard X-ray/ soft gamma- 
ray regime [9£]. The GeV spectrum during flares in- 
dicates that the emission is confined to the synchrotron 
component of the SED, a conclusion which is supported 
by the rapid timescale of the events (since the inverse 
Compton or Bremsstrahlung cooling time of the emit- 
ting electrons is much greater than the observed flare 
duration). At higher energies, some evidence for an en- 
hanced flux during HE flare states has been presented by 



ARGO-YBJ ilOOll . I ACT measurements do not support 
these results, but are not necessarily in conflict, given 
the differing duty cycles. Detailed measurements with 
lACTs during future flare states are required to resolve 
this question. 

5.3. Pulsar Wind Nebulae 

Pulsar wind nebulae are the most abundant class of 
known VHE emitters in the Galaxy, with ~ 30 firm ex- 
amples, and numerous other sources where the PWN 
association is more tentative (for reviews see e.g. flOll 
102ill03.1 ). The essential emission mechanisms - shock 
accelerated leptons producing synchrotron and inverse 
Compton radiation - have already been described for the 
case of the Crab PWN, but the Crab is far from the typi- 
cal object. Understanding of the structure and evolution 
of PWN has advanced significantly over the past few 
years, in particular thanks to the high resolution X-ray 
imag ing provided by Chandra (see Gaensler and Slane 
lll04! l for a detailed review). Initially, the PWN expands 
uniformly from the central pulsar, while at later stages 
the nebula may be confined and distorted by the reverse 
shock from the expanding supernova remnant (SNR). 

At TeV energies, young PWN are usually still embed- 
ded within their parent SNR and are point-like, within 
the angular resolution of lACTs. They are position- 
ally coincident with a bright X-ray synchrotron nebula 
powered by a pulsar w ith v ery high spin-down lumi- 
nosity (e.g. G0.9H-0.1 Col, HESS J1813-178 iflOfill . 
G54.1H-0.3 HotIi '). More evolved PWN, with ages > 
10,000 years, are usually much larger, and their TeV 
emission can be spatially resolved and mapped. The 
pulsar powering the nebula is often offset from the cen- 
ter of the TeV emission, probably for reasons related to 
densi ty gradients in the medium surrounding the SNR 
1 10811 . Remarkably, the TeV nebulae are often two or 
three orders of magnitude larger than the corresponding 
X-ray PWN, and the TeV PWN sizes tend to increase 
with age, while the X-ray PWN sizes show the opposite 
trend. This can be understood as a result of the fact that 
the electron population which is responsible for the TeV 
inverse Compton flux has lower energies than the elec- 
trons which produce the X-ray synchrotron emission. 
They therefore cool more slowly, and survive for longer 
so, while the X-ray nebula is dominated by freshly ac- 
celerated particles, the TeV nebula can record the entire 
history of particle propagation away from the pulsar. 
A natural result of this is that the observed TeV spec- 
trum should vary with distance from the pulsar One of 
the best examples of this is HESS J 1825- 137, associ- 
ated with the PWN of the pulsar PSRJ 1826- 1334 iToj. 
Figure [12] shows the spatially dependent spectra for this 
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Figure 1 1 : Left: The spectral energy distribution of the Crab Nebula from Abdo et al. (90ll . Their fit to the synclirotron component is shown (blue 
dashed line), as well as inverse Compton spectra from Atoyan and Aharonian |97] for assumed magnetic field strengths of 100 /iG (solid red line), 
200 /jG (dashed green line) and 300 fjG (dotted blue line). Right: High Energy spectrum of the Crab Pulsar The black dashed line shows a fit of 
a power law with exponential cut-oft' to the Fermi-hAT data alone; the solid fine shows a broken power law fit to the LAT and lACT data. The red 
dashed line and bowtie shows a power law fit to the VERITAS points alone. Figure courtesy of N. Otte (priv. comm.). 



source, which soften with increasing distance from the 
pulsar This is interpreted as the natural effect of both 
inverse Compton and synchrotron cooling of the elec- 
tron population during propagati on. A counter-example 
is the case of Vela-X UM I 

111 in . in which no spectral 
variability is seen over the extended nebula, suggesting 
that cooling does not play an important role. 

5.4. Supernova Remnants 

The search for the origin of the cosmic rays trig- 
gered the development of gamma-ray astronomy, and 
continues to motivate many gamma-ray observations. 
Chief among these is the study of supernova remnants, 
which are believed to efficiently accelerate particles at 
the shock front where the expanding SNR encounters 
the surrounding medium (e.g. ill 12ll ). This hkely occurs 
through diffusive shock acceleration (first order Fermi 
acceleration), in which charged particles are reflected 
from magnetic inhomogeneities and repeatedly cross 
the s hock fron t, gaining energy with each crossing (see 
e.g. ll 13lll 1411 ). As well as plausibly providing enough 
energy to explain the observed Galactic cosmic ray pop- 
ulation, this process naturally produces a power law dis- 
tribution of particle energies with an index of ~ 2, which 
matches the cosmic ray spectrum (after accounting for 
diffusion and escape). In recent years, the importance 
of magnetic field amplification by the accelerated par- 
ticles themselves has been increasingly recognised, and 
plays a particular role in explaining the existence of the 
highest energy Galactic cosmic rays, around the cosmic 
ray knee region (at ~ 3 x 10'^ eV). 

The evidence for efficient leptonic acceleration in 
SNRs is now clearly established (e.g. [115]); how- 
ever, the question of whether SNR are efficient hadron 



accelerators is more difficult to answer A definitive 
measurement would be the detection of high energy 
neutrinos from an SNR, but the expected fluxes are 
likely below the sensitivity thresholds of current neu- 
trino observatories. Gamma-ray observations may pro- 
vide the key, since the interactions of high energy nu- 
clei with target material produce neutral pions, which 
decay immediately into gamma-rays. Disentangling the 
spectral signature of this process from other sources of 
gamma-ray emission (i.e. leptonic inverse Compton and 
bremsstrahlung processes) is difficult, but not impossi- 
ble. Two classes of gamma-ray source are of interest 
for these studies: those which can be clearly associated 
with SNR shells, based on the gamma-ray morphology, 
and sources which are coincident with a massive volume 
of target material, such as molecular clouds. 

The definitive association of gamma-ray emission 
with an SNR shell is often difficult to make, due to 
the presence of other potential counterparts, particu- 
larly PWN. A handful of shell- type SNRs have been un- 
equivocally identified as gamma-ray sources by lACTs. 
This identification can be made on the basis of the 



observed sheU morphology, (RXJ 1713.7-3946 Illl6ll . 
RXJ 0852.0-4622 (Vela Jr) jjUl, HESS J173 1-3467 
ifTTs fl. SN1006 iim and, possibly, RCW 86 fl2^), or, 
in the case of Tycho's SNR [ 12 IJ, on the positional co- 
incidence, coupled with the fact that the progenitor was 
a known Type la exposion, and so no compact object 
is present. Figure [13] shows the gamma-ray map for 
flie first SNR shefl to be resolved, RXJ 1713.7-3946. 
The recent addition of Ferm/-LAT observations to the 
broadband spectrum of RXJ 1713.7-3946 [122] are con- 
sistent with a leptonic origin as the dominant mecha- 
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Figure 12: Inset: H.E.S.S. gamma-ray excess map for HESS J1825- 
137. The wedges show the radial regions with radii in steps of 0.1° 
in which the energy spectra were determined. The main figure shows 
the differential energy spectra for the regions illustrated in the inset, 
scaled by powers of 10 for the purpose of viewing. The spectrum for 
the analysis at the pulsar position is shown as a re ference along with 
the other spectra as dashed line. Figure from |109|| : see that paper for 
full details. 



nism for the gamma-ray emission. A counter-example 
is illustrated by the spectrum in Figure [13] which cor- 
responds to Tycho's SNR. In this case, the fact that the 
broadband gamma-ray spectrum can be fit with a sin- 
gle hard power la w from 50 MeV to 10 TeV favours 
a hadronic origin 1 1231 124 1 . This interpretation is not 
completely compelling, however, given the large statis- 
tical errors in the measurements, and the impact of vari- 
ous unknown parameters such as the SNR distance, and 
possible enhancements of the gamma-ray flux due to a 
nearby molecular cloud. Additionally, Atoyan and Der- 
mer 112511 describe a two-zone leptonic model which 
provides an acceptable spectral fit. Future measure- 
ments of the spectrum below 500 MeV, and deeper ex- 
posure at TeV energies, will further test the differing 
interpretations. 

The intensity of gamma-ray emission due to hadronic 



interactions depends upon the flux of high energy nu- 
clei, and also upon the density of target material. Re- 
gions of high matter density (e.g. molecular clouds with 
densities > 100 cm ■^), situated close to sites of particle 
acceleration (such as SNRs), can therefore be expected 
to produce a large gamma-ray flux due to hadronic in- 
teractions. At TeV energies, ~ 10 likely candidates for 
this process have been identified. The task of identifica- 
tion is complicated, both for the usual reasons of source 
confusion, and also because the evidence for a molec- 
ular cloud / SNR interaction is only definitive in those 
cases where the cloud morphology is visibly deformed 
by the expanding SNR, and/or where sites of hydroxyl 
(OH) maser emission indicate the presence of shocked 
molecular material. One of the best examples of this 
source class is the old remnant W28 lll26ll . H.E.S.S. 
observations of this region show four distinct sites of 
emission, with three of the four showing a resolvable 
angular extent (~ 10')- Each of the TeV sources is po- 
sitionally coincident with a molecular cloud. Assuming 
that the TeV emission is due to hadronic cosmic rays in- 
teracting with the cloud material, the cosmic ray density 
is inferred to be a factor of 10 to 30 times greater than 
in the solar neighbourhood. 

5.5. Star forming regions 

The process of diff'usive shock acceleration is not lim- 
ited to supernova remnant shells. An alternative sce- 
nario invokes particle acceleration at the shock formed 
by the collision between the supersonic stellar winds 
of massive stars in close binary systems. Stellar winds 
may also become collectively important in large assem- 
blies of massive stars. The combined effect of the stellar 
winds, coupled with the effect of multiple SNRs, results 
in an overall wind from the cluster which forms a giant 
superbubble in the interstellar medium. Particle accel- 
eration can occur where the cl uster wind interacts with 
the surrounding medium (e.g. 

Massive star associations are naturally likely to host 
other potential source counterparts for TeV emission, 
such as compact object binary systems, individual su- 
pernova remnants and pulsar wind nebulae. A case in 
point is the young, open cluster Westerlund 2, contain- 
ing the Wolf-Rayet binary system WR 20a. This was 
originally suggested as a plausible counterpart to the 
unidentified source HESS J1023-575, with the emis- 
sion presumed to be connected to either t he W olf-Rayet 
binary, or the combined cluster wind lll28ll . A re- 
assesment of this region, informed by a deeper H.E.S.S. 
exposure and results from ferm/-LAT, alters the picture 
[129]. The LAT detects an energetic pulsar, PSR J1022- 
5746, which drives a PWN which is bright in GeV 
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Figure 13: Left: H.E.S.S. map of gamma-ray excess eve nts fo r RXJ 1713.7-3946 - the first SNR shell to be resolved at TeV energies. Figure from 
llltf l. Right: The broadband SED of Tycho's SNR from |124|| . together with models for the various emission components (dominated by hadronic 
processes in the gamma-ray band). See paper for details. 



gamma-rays il30ll . Given the ubiquity of bright TeV 
PWN, this now seems the most likely explanation for 
the TeV source. A similar conclusion may arise for 
the first TeV source to be linked with a massive star 
association, TeV 2032-1-4130 (coincident with the Cyg 
OB2 association). In this case, the LAT pulsar (PSR 
J2032-I-4127) is sufficiently energetic to explain the TeV 
emission, although no PWN has been detected as yet. 

Other unidentified sources which have been linked 
with massive star clusters and associations include 
HESS J1646-458 (Westerlund 1) jlJ^, HESS J1614- 
518 (Pismis 22) 1132], HESS J1848-018 (W43, which 
hosts Wolf-Rayet starWR121a) [133] and W49A ifBl. 
a massive star forming region. For all of these, however, 
the evidence that particle acceleration in stellar winds is 
the driving force behind the gamma-ray emission is not 
definitive (e.g. lil35ll '). 

5.6. Compact Object Binary Systems 

Despite many early unconfirmed claims, the first 
definitive detection of a TeV gamma-ray binary system 
was not published until 2005. The population has grown 
slowly, and now consists of four clearly identified sys- 
tems, plus marginal evidence for transient emission as- 
sociated with Cyg X-1 The gamma-ray emission 
from binaries is believed to be powered either by accre- 
tion (most likely onto a black hole), or by a pulsar wind. 
In the case of accretion, partic le acceleration takes place 
in relativistic jets (e.g. lll37ll ). In the pulsar wind sce- 
nario, the acceleration occurs either in shocks where the 
pulsar wind encounters the circumstellar environment 



would be extremely important, since there is no doubt 
that this system hosts a black hole. This is in contrast 
to all of the other TeV binaries, in which the compact 
object is either known to be, or may be, a neutron star 
Here we briefly summarize the results for each of the 
four well-studied objects. 

PSR B1259-63/LS 2883: This was the first gamma- 
ray binary system to be firmly detected at TeV energies, 
and the first known variable VHE source in our Galaxy 
lll4Ql . The system comprises a 48 ms pulsar orbiting a 
massive B2Ve companion. The orbit is highly eccentric 
(e = 0.87), with a period of 3.4 years. The TeV emis- 
sion exhibits two peaks, approximately 15 days before 
and after periastron. Various authors have attempted 
to explain the double bumped VHE lightcurve within 
a 'hadronic disk scenario', in which a circumstellar disk 
provides target material for accelerated hadrons, lead- 
ing to 7T^ productio n and subsequent TeV gamma-ray 



emission 111411 114211 . The 2007 H.E.S.S. observations 



(e.g. 
itself 



138] ), or possibly within the pulsar wind zone 
13911 . The detection of Cyg X-1, if confirmed. 



disfavour this, since the onset of TeV emission occurs 
~ 50 days prior to periastron, well before interactions 
with the disk could be expected to play a significant role 
1 14311 . Leptonic scen arios have also been discussed in 
e.g. Kangulyan et al. 11 14411 . 

The recent discovery of an extended and variable 
radio structure in PSR B1259-63/LS 2883 at phases 
far from periastron provides definitive evidence that 
non-accreting pulsars orbiting massive stars can pro- 
duce variable and extended radio emission at AU scales 
1 14511 ■ This is important, since similar structures in 
LS 5039 and LS I -h6r303, where the nature of the 
compact object is not certain, have been used to argue 
for the existence of jets driven by accretion onto a black 
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Figure 14: The H.E.S.S. flux (bottom) and photon index (top) for 
LS 5039 as a function of orbital phase. Figure from Aharonian et 
al. fT4l. 



hole. 

LS 5039: LS 5039 consists of a compact object, either 
neutron star or black hole, orbiting a massive 06.5V 
(~23 Mq) star in a 3.9 day orbit. Observations by 
H.E.S.S. in 2004 revealed that LS 5039 is a bright 
source of VHE gamma-rays lfl46t . Unlike PSR B1259- 
63 (and, to a lesser extent, LS I +61°303) LS 5039 is al- 
most perfectly suited to TeV observations, with a short 
orbital period and a convenient declination angle, al- 
lowing sensitive observations at all phases over numer- 
ous orbits. The VHE emission measured by H.E.S.S. 
is modulated at the orbital period, peaking around in- 
ferior conjunction, when the compact object is clos- 
est to us and co-aligned with our line-of-sight (Fig- 
ure [T4]i. The spectrum is also orbitally modulated, ap- 
pearing significantly harder around inferior conjuction 
(F - L85 + 0.06stat ± O.lsyst), but with an exponential 
cut-off at Eo = 8.7 + 2.0 TeV. At GeV energies, the 
source is detected by Fenni-LAI at all orbital phases, 
with the emission peaking close to superior c onjuc tion, 
in apparent anti -phase with the VHE results lll47ll . A 
sharp spectral cut-off at Eg - 1.9 + 0.5 GeV is observed 
in the LAT data near superior conjunction, indicating 
that the VHE spectra cannot be simply a smooth extrap- 
olation of the lower energy emission. 

LS I +61°303: Similar- to LS 5039, LS I +61°303 con- 
sists of a compact object, either neutron star or black 
hole, in this case orbiting a BOVe star with a circum- 
stellar disk (~12.5 Mq) in a 26.5 day orbit. The de- 



tection of a variable VHE source at the location of 
LS I +61°303 with MAGIC fu^, later confirmed by 
VERITAS 1.1491 . established this source as a gamma- 
ray binary. The object is now one of the most heav- 
ily observed locations in the VHE sky, with deep ex- 
posures by the two observatories spread over half a 
decade. Despite this, the VHE source is much less well- 
characterized than LS 5039, owing to its relatively weak 
VHE flux, and an inconvenient orbital period which 
closely matches the lunar cycle, making observations 
over all orbital phases almost impossible within a sin- 
gle observing season. VHE emission was originally de- 
tected close to apastron, between phases <p = 0.5 - 0.8. 
In contrast to this, VHE observations between 2008 and 
2010 showed that, at least during the orbits when the 
source was observed, the apastron fl ux was mu ch lower 
than during the previous detections [150i 151 1. The de- 
tection of TeV emission by VERITAS during a single 
episode close to superior conjunction complicates the 
picture even further As with LS 5039, the ferm/-LAT 
GeV emission peaks closer to periastron, and the spec- 
trum displays sharp cut-off, at Eg = 6.3 + LI GeV. Long 
term variability in the GeV band has also been observed 

iliii. 

HESS J0632+057: This TeV source was serendip- 
itously detected during H.E.S.S. observations of the 
Monoceros Loop SNR region, and noted as a poten- 
tial b inary primarily because of its small angular extent 
lll53ll . Subsequent observations revealed it to be a vari- 
able TeV source 1 154.1 . and co-located with a variable 
radio and X-ray sourc e, at the position of a massive Be 
star, MWC 148 fl5?l. X-ray observations with Swift 
recently provided definitive proof of its binary nature, 
with the measurement of a 321 + 5 day periodicity in the 
lightcurve |156]. The TeV light curve displays a broad 
gamma-ray flare close to th e X-ray maximum, with a 
duration of ~ 40 days lll57ll . No GeV source has been 
detected. 

A number of competing processes likely contribute 
to the variability observed in the gamma-ray binaries. 
In particular, the efficiency of inverse Compton gamma- 
ray production, as well that of VHE gamma-ray absorp- 
tion (through pair production), changes as a function 
of orbital phase. This can go some way towards ex- 
plaining the apparent phase shift between the GeV and 
TeV lightcurves, for example in the case of LS 5039. 
There are clearly other effects which contribute, how- 
ever, as demonstrated by the long term instability of 
LS I +61°303. The sharp GeV cut-offs in these sys- 
tems are also difficult to explain, and may indicate that 
the GeV and TeV emission components do not have the 
same origin. 
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A final comment should be reserved for the Fermi- 
LAT source 1 FGL J1018.6-5856, which was recently 
identified as a Ge V bi nary system, with a period of 
16.58 + 0.04 days lll58ll . This source resides in a com- 
plex region, coincident with the center of the SNR 
G284.3-1.8, and close to a LAT pulsar (at a distance 
of 35'). The extended, unidentified HESS source, 
HESS J1018-589, overlaps with the GeV binary loca- 
tion, and appears to consi st of a point-like source over- 
laid on a diffuse structure lll59ll . While this is indicative 
of a new TeV binary there is, as yet, no evidence for 
variability in the point-like emission, making the identi- 
fication still uncertain. 

5. 7. Globular Clusters 

A single globular cluster, Terzan 5, has been sug- 
gested as the prob able counterpart of a TeV source 
(HESS J1747-248) 1160]. If the association is correct, 
the emission is likely related to the large population of 
millisecond pulsars in this cluster, wh ich p rovide an 
injection source of relativistic leptons [161]. Inverse 
Compton gamma-ray photons result when these elec- 
trons upscatter low energy photons of the intense stellar 
radiation field. The H.E.S.S. source is extended, and 
slightly offset from the cluster core (although there is 
significant overlap). The probability of this being sim- 
ply a chance positional coincidence is ~ 10 

Globular clusters have also been favored targets in 
searches for dark matter particle annihilation signa- 
tures, since they may have been generated in dark mat- 
ter mini-haloes before the formation of galaxies took 
place , and thus retain a significant dark matter compo- 
nent f 162 , 163 1 . Limits have been placed on NGC (5388 , 
M15, Omega Centauri, 47 Tuc, M13, and M5 ( lfl64ll . 
and references therein). 

5.8. Unidentified objects 

lACTs are able to locate point sources with reason- 
ably good accuracy (typically < 1' for a moderately 
strong source). Extragalactic TeV sources can therefore 
usually be firmly identified with a single counterpart at 
other wavelengths. For Galactic sources, identification 
poses more of a problem, and around one third of the 
Galactic sources lack a firm identification. While the 
diffuse Galactic background emission, which dominates 
at GeV energies, is not significant, the TeV sources 
themselves are mostly extended, and can often be plau- 
sibly associated with multiple counterparts. PWN, in 
particular, often have their brightest TeV emission off- 
set from the location of the parent pulsar or X-ray PWN, 
which may not yet have been detected (a number of new 



energetic pulsars have been located in follow-up obser- 
vations of unidentified TeV sources). In other cases, de- 
spite deep X-ray and radio follow-up observations, no 
reasonable counterpart has yet been found. 

The unidentified Milagro sources also pose some 
interesting questions. While lACTs have identified 
sources associated with these objects, the TeV sources 
are typically much smaller in angular extent, and can- 
not account for the entire Milagro flux. One possible 
explanation for this is that there is a diffuse high energy 
component to the emission, which is difficult to resolve 
with lACTs, given their relatively small fields of view. 

6. The Future 

The field of ground-based gamma-ray astronomy has 
expanded dramatically over the past 10 years, but it is 
worth noting that the observatories currently operating 
are far from reaching the physical limits of the detec- 
tion techniques. We therefore conclude this review with 
a brief discussion of some of the instrumental develop- 
ments which can be expected over the next decade. 

The lACTs currently operating have all made signifi- 
cant efforts to maintain, and improve, sensitivity since 
they were first commissioned. H.E.S.S. recoated the 
telescope mirrors, and successfully developed sophis- 
ticated analysis tools with which to exploit the data. 
VERITAS relocated their original prototype telescope 
to provide a more favorable array layout, halving the 
time required to detect a weak source. Most signifi- 
cantly, MAGIC added a second telescope of similar de- 
sign to the first, providing a stereo pair with a base- 
line of 85 m. All of these observatories have further 
upgrade plans. Both MAGIC and VERITAS are imple- 
menting camera upgrades - in the case of VERITAS this 
involves the replacement of all of the camera PMTs with 
more sensitive, super-bialkali devices in summer 2012. 
H.E.S.S. are constructing H.E.S.S. II - the addition of a 
single large telescope, with 600 m^ mirror area, to the 
centre of the array. Figure [T5]shows the telescope struc- 
ture in November 2011 . The mirrors and camera will be 
installed in the first half of 2012. 

Novel approaches are also being pursued by many 
smaller projects. A particularly nice example of this is 
FACT (the First G-APD Cherenkov Telescope), which 
has recently demonstrated the application of Geiger- 
mode Avalanche Photodiodes for Cherenkov astron- 
omy. The FACT telescope consists of a 1440-pixel G- 
APD camera at the focus of one of the original HEGRA 
telescopes. Geiger-APDs hold great promise as a poten- 
tial replacement for PMTs, since they are robust, and of- 
fer much better photon conversion efficiency. Figure [16] 
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Figure 15: The steel structure for the 600 H.E.S.S. II telescope. 
MiiTors and camera will be installed in the first half of 2012. (Note 
that the H.E.S.S. II telescope is located in between the two H.E.S.S. I 
telescopes shown in the image, not in the foreground). Figure courtesy 
of the H.E.S.S. collaboration http://www.mpi-hd.mpg.de/hfm/HESS/ 



shows some "first light" images from FACT. Other 
projects in development include GAW (Gamma Air 
Watch), MACE (Major Atmospheric Cerenkov Tele- 
scope Experiment) and LHASSO (Large High Altitude 
Air Shower Observatory). LHASSO is an ambitious 
project which will be located near the site of the ARGO- 
YBJ experiment in Tibet. It is planned to consist of four 
water Cherenkov detectors, two lACTs, three fluores- 
cence telescopes and a large scintillator array. 

For the particle detectors, the next stage in instrumen- 
tation is HAWC (the High Altitude Water Cherenkov 
Observatory). HAWC will consist of 300 individual wa- 
ter Cherenkov tanks, at an altitude of 4100 m in Mexico. 
The final array is expected to be 15 times more sensitive 
than Milagro, and will be a powerful tool for surveying, 
and for observations of transient phenomena. A proto- 
type system is already in operation, and science opera- 
tions will start in spring 2012, with completion of the 
full array expected in 2014. 

The most ambitious future project is CTA (the 
Cherenkov Telescope Array). This is described in de- 
tail in 11651]. Briefly, it comprises an array of imaging 
atmospheric Cherenkov telescopes covering ~ 1 km-, 
providing a factor of 5-10 improvement in sensitivity 
in the lOOGeV-lOTeV range, and extending the en- 
ergy range both above and below these values. Both 
a northern and a southern site are envisaged, and the ar- 



Figure 16: First light cosmic ray images from FACT (First G-APD 
Cherenkov Telescope). Four different events are shown; the upper 
right event shows the characteristic ring image produced by a local 
muon. The camera contains 1440 Geiger-mode avalanche photodi- 
odes, installed on one of the original HEGRA telescopes at the Roque 
de los Muchachos on La Palma. Figure courtesy of the FACT collab- 
oration http://fact.ethz.ch/first 



ray will be operated as an open observatory. Multiple 
telescope designs are planned, including small (few m), 
medium (10 - 15 m) and large (20 - 30 m) diameter re- 
flectors, as well as two-mirr or tel escope designs, such 
as the Schwarzchild-Couder lll66ll . 

In conclusion, TeV gamma-ray astronomy now de- 
scribes a broad astronomical discipline which addresses 
a wide, and expanding, range of astrophysical topics. 
With planned instrumental developments, it is not un- 
reasonable to expect the source catalogue to exceed 
1000 objects within the next decade. Much of the TeV 
sky remains relatively unexplored. Less than ~ 10% of 
the sky has been observed with ~ 10 milUCrab sensi- 
tivity at 1 TeV, and the sensitive exposure to transient 
events and widely extended sources at these energies is 
much lower The likelihood of continued exciting re- 
sults is certain, both for the known sources and source 
classes, and for new discoveries. 
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